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A mechanism of disappearance and formation of the Efimov levels of 
the helium 4 He3 trimer is studied when the force of interatomic interaction is 
changed. The resonances including virtual levels are calculated by the method 
based on the solving the boundary value problem, at complex energies, for 
the Faddeev differential equations describing the (2 + 1— ►2 + 1; 1 + 1 + 1) 
scattering processes. 

I. INTRODUCTION 

The system of three 4 He atoms is of considerable interest for various fields of physical 
chemistry and molecular physics. The present paper is a sequel of studies of the 4 He3 
system undertaken in jlj, where the repulsive component of the He-He interaction at short 
distances between atoms is approximated by a hard core. This allows one to investigate 
the 4 He3 system within a mathematically rigorous method of solving a three-body problem 
in the Boundary-Condition Model developed in 0|§. In 0, such an approach has been 
successfully applied for calculating not only scattering but also the ground- and excited- 
state energies of the helium trimer. Investigation made in [[I] has shown that the method 
proposed in 0|| is well suited for performing three-body molecular computations in the 
case where repulsive components of interatomic interactions are of a hard core nature. 

There is a series of works |4],|,[l| showing that the excited state of the 4 He trimer is 
initiated indeed by the Efimov effect ||. In these works the various versions of the Aziz 
4 He- 4 He potential were employed However, the basic result of Refs. on the excited 

state of the helium trimer is the same: this state disappears when the interatomic potential 
is multiplied by the "amplification factor" A of order 1.2. More precisely, if this potential 
is multiplied by the increasing factor A > 1 then the following effect is observed. First, 
the difference e<j(A) — e[ 1 \\) between the dimer energy e d (X) and the energy of the trimer 
excited state ^"^(A) increases. Then the behavior of this difference radically changes and 
with further increase of A it monotonously decreases. At A ~ 1.2 the level E^ disappears. 
It is just such a nonstandard behavior of the energy E\ [ (A) as the coupling between helium 



* Contribution to Proceedings of the First International Conference on Modern Trends in Compu- 
tational Physics, June 15-20, 1998, Dubna (Russia). LANL E-print |physics/9810005| . 

^On leave of absence from the Laboratory of Theoretical Physics, Joint Institute for Nuclear 
Research, Dubna, 141980, Russia 



1 



atoms becomes more and more strengthening, points to the Efimov nature of the trimer 
excited state. And vice versa, when A slightly decreases (no more than 2%), the second 
excited state E^f 1 appears in the trimer @,|]. 

Here we present the results of our numerical study of a mechanism of disappearance 
and formation of the Efimov levels of the helium 4 He3 trimer using the method of search 
for resonances in a three-body system on the basis of the Faddeev differential equations. 
The idea of the method formulated and proved in |7j consists in calculating the analytic 
continuation of the component So(z) of the scattering matrix corresponding to the (2 + 1 — > 
2 + 1) process in the physical sheet. For the potentials we use, the three-body resonances 
(including virtual levels) lying in the unphysical sheet of energy z plane adjoining the physical 
sheet along the interval (e d , 0) are the roots of the function S (z) in the physical sheet. We 
have earlier employed this method for computing resonances as roots of So {z) in the three- 
nucleon problem ||. 



II. METHOD 

In this work we consider the three-atomic 4 He system with the total angular momentum 
L = 0. The angular partial analysis reduces the initial Faddeev equation for three identical 
bosons to a system of coupled two-dimensional integro-differential equations (see Ref. |]| 
and references therein) 



d 2 d 2 ( 1 1 \ 



Fi(x,y) (1) 



dx 2 dy 2 \x 2 y 2 

= { -V(x)^i(x,y), x > c 
~ [ 0, x < c. 

Here, x, y stand for the standard Jacobi variables and c, for the core range. At L = 
the partial angular momentum I corresponds both to the dimer subsystem and a com- 
plementary atom. The energy z can get both real and complex values. The He-He po- 
tential V(x) acting outside the core domain is assumed to be central. The partial wave 
function *$?i(x, y) is related to the Faddeev components Fi(x,y) by ^i(x,y) = Fi(x,y) + 
£i' f^dr)h w (x,y, rf) F v {x' , y') where x' = (\x 2 + § y 2 - & xyr]) 1 / 2 ,y' = (\x 2 + \ y 2 + 
^ xyr]) 1 / 2 and 1 < rj < 1. The explicit form of the function hui can be found in Refs. p|,|T0 . 
The functions F[(x,y) satisfy the boundary conditions 

Fi(x, y) \ x=0 = Fi(x, y) \ y=0 = and %(c, y) = . (2) 

Here we only deal with a finite number of equations ([!]), assuming that I < Z max where Z max 
is a certain fixed even number. The condition < / < / max is equivalent to the supposition 
that the potential V(x) only acts in the two-body states with I = 0,2, ... , Z max . We assume 
that the potential V(x) is finite, i. e., V(x) = for x > r , r > 0. The asymptotic conditions 
as p — > oo and/or y — » oo for the partial Faddeev components of the (2 + 1 — > 2 + 1 ; 1 + 1 + 1) 
scattering wave functions for z = E + iO, E > 0, read (see, e. g., Ref. 0) 

F t (x,y;z) = 5i i/; d (x) {sm(^z - e d y) + exp(iV^ - t d y) [& (z) +o(l)]} 

+ ex E M^ M)+o(1)] . (3) 
VP 
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We assume that the 4 He dimer has an only bound state with an energy e^, < 0, and 
wave function ^(x), i>d{ x ) — for < x < c. The notations p, p = ^/x 1 + y 2 , and 
6,6 = arctan(?//a;), are used for the hyperradius and hyperangle. The coefficient ao(z), 
z = E + iO, for E > is the elastic scattering amplitude. The functions Ai(E + i0,6) 
provide us, at E > 0, the corresponding partial Faddeev breakup amplitudes. For real 
z = E + iO, E > €d, the (2 + 1— >-2 + 1) component of the s-wave partial scattering matrix 
for a system of three helium atoms is given by the expression 

S o (z) = l + 2iao02). 

Our goal is to study the analytic continuation of the function Sq(z) into the physical sheet. 
As it follows from the results of Refs. 0, the Sq(z) is just that truncation of the total 
scattering matrix whose roots in the physical sheet correspond to location of resonances in 
the unphysical sheet adjoining the physical one along the spectral interval (e^, 0). 



III. RESULTS OF COMPUTATIONS 

In the present work we make use of the Faddeev equations ([!]) considered together with 
the boundary conditions (fj), @ to calculate the values of the 4 He 3 scattering matrix S (z) 
in the physical sheet. We search for the resonances including the virtual levels as roots of 
Sq(z) and for the bound-state energies as positions of poles of S (^). All the results presented 
below are obtained for the case Z max = 0. 

In all our calculations, h 2 /m = 12.12 KA 2 . As the interatomic He - He - interaction we 



employed the widely used semiempirical potential HFD-B |L2] . The value of the core range c 
is chosen to be equal 1.3 A providing at least six reliable figures of the dimer binding energy 
€d and three figures of the trimer ground state energy i?f . A detailed description of the 
numerical method we use is presented in Ref. JT|. In contrast to [|lj, in the present work we 
solve the block-three-diagonal algebraic system on the basis of the matrix sweep method. 
This allows us to dispense with writing the system matrix on the hard drive and to carry out 
all the operations related to its inversion immediately in RAM. Besides, the matrix sweep 
method reduces almost by one order the computer time required for computations on the 
grids of the same dimensions as in jl). 

We searched for the resonances (roots of the function Sq(z) on the physical sheet) and 
bound-state energies (roots of the function Sg (z) for real z < e<j) of the helium trimer by us- 
ing the complex version of the secant method. We found positions of the four "resonances", 
the roots of S (^), in case of the grid parameters N e = N p = 600 and p max = 600 A. Complex 
roots of the function S (z) are located at points (—2.34 + i0.96)mK, (—0.59 + i2.67)mK, 
(2.51 + i 4.34) mK and (6.92 + i 6.10) mK. These "resonances" are situated beyond the do- 
main of scattering matrix holomorphy where the applicability of our method is proved 
Till - So we do not consider the roots of function Sq(z) as genuine resonances for the 4 He3 



system. However it is remarkable that the "true" (i. e., getting inside U^) virtual levels and 
then the energies of the excited (Efimov) states appear just due to these (quasi) resonances 
when the potential V(x) is weakened. 

Following mHI]], instead of the initial potential V(x) = Vhfd-b(^), we consider the po- 
tentials V(x) = X- Vhfd-b^)- To establish the mechanism of formation of new excited states 
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in the 4 He trimer, we have first calculated the scattering matrix So (z) for A < 1. We have 
found that for a value of A slightly smaller than 0.9885, the (quasi)resonance closest to the 
real axis gets on it and transforms into a virtual level of the second order. This virtual level 
is preceded by the (quasi) resonances z = (—1.04 + i 0.11) mK (^/|e^| = —1.58 + i 0.168) for 
A = 0.989 and z = (-0.99 + i 0.04) mK (z/\e d \ = -1.59 + i 0.064) for A = 0.9885. The origi- 
nating virtual level is of the second order since simultaneously with the root of the function 
So(z), a l so the conjugate root of this function gets on the real axis. With a subsequent 
decrease of A the virtual level of the second order splits into a pair of the virtual levels of 
the first order which move in opposite directions. One of the virtual level moves towards 
the threshold and "collides" with it at A < 0.98. For A = 0.975 the function S (z) instead 
of the root acquires a new pole corresponding to the second excited state of the trimer with 
the energy E[ . 

To determine the mechanism of disappearance of the excited state of the helium trimer 
we have calculated the scattering matrix So(z) when the two-body interactions become 
stronger owing to the increasing coupling constant A > 1. We found that this disappearance 
proceeds just according to the scheme of the formation of new excited states; only the order 
of occurring events is inverse. In the interval between A = 1.17 and A = 1.18 there occurs 
a "jump" of the level on the unphysical sheet and it transforms from the pole of the 
function S (z) into its root corresponding to the trimer virtual level. 



ACKNOWLEDGEMENTS 

The authors are grateful to Prof. V. B. Belyaev and Prof. H. Toki for help and assistance 
in calculations at the supercomputer of the Research Center for Nuclear Physics of the Os- 
aka University, Japan. One of the authors (A. K. M.) is much indebted to Prof. W. Sandhas 
for his hospitality at the Universitat Bonn. The support of this work by the Deutsche 
Forschungsgemeinschaft and Russian Foundation for Basic Research is gratefully acknowl- 
edged. 



[1] E. A. Kolganova, A. K. Motovilov, S. A. Sofianos: J. Phys. B. 31 (1998) 1279 (LANL E-print 

|physics/9612012| ). 
[2] S. P. Merkuriev, A.K. Motovilov: Lett. Math. Phys. 7 (1983) 497. 

[3] S. P. Merkuriev, A. K. Motovilov, S. L. Yakovlev: Theor. Math. Phys. 94 (1993) 306 (also see 

LANL E-print |iucl-th/9606022| ). 
[4] B. D. Esry, C. D. Lin, C. H. Greene: Phys. Rev. A. 54 (1996) 394. 
[5] T. Cornelius, W. Glockle: J. Chem. Phys. 85 (1986) 3906. 
[6] V. Efimov: Nucl. Phys. A. 210 (1973) 157. 

[7] A. K. Motovilov: Math. Nachr. 187 (1997) 147 (LANL E-print [funct-an/9509003l ). 
[8] E. A. Kolganova, A. K. Motovilov: Phys. Atom. Nucl. 60 (1997) 177 (LANL E-print 
|rmcl-th/9602001| ) ; also see LANL E-print |nucl-th/9702037 . 



4 



[9] L. D. Faddeev, S. P. Merkuriev: Quantum Scattering Theory for Several Particle Systems. 
Doderecht: Kluwer Academic Publishers 1993. 
[10] S. P. Merkuriev, C. Gignoux, A. Laverne: Ann. Phys. (N.Y.) 99 (1976) 30. 
[11] E. A. Kolganova, A. K. Motovilov: Preprint JINR E4-98-243 (LANL E-print 
|physics/9808027| ). 

[12] R. A. Aziz, F. R. W. McCourt, C. C. K. Wong: Mol. Phys. 61 (1987) 1487. 



5 



